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(54) Multi-rate digital filter for audio sample-rate conversion 



(57) A method for providing a sample-rate conver- 
sion ("SRC") filter on an input stream of sampled data 
provided at a first rate, to produce an output stream of 
data at a second rate different from the first rate. The 
input stream of sampled data is operated on with a first 
low-order interpolation filter routine 42 to produce a first 
stream of intermediate data The first stream of interme- 
diate data is operated on with a first simplified interpola- 



tion filter routine 44,46, having a substantially small 
number of operations to calculate the coefficients 
thereof, to produce a second stream of intermediate 
data. The second stream of intermediate data is oper- 
ated on with a first decimating filter routine 48 to pro- 
duce the output stream of data 
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Description 

Technical Field of the Invention 

5 [0001] This invention relates to audio sample-rate conversion systems, and more particularly relates to multistage 
sample-rate conversion filters. 

Background of the Invention 

10 [0002] For both historical and technical reasons, there have existed a number of industrial standards on audio dig- 
ital signal sample rates. The well known examples are the 44.1 kHz sample rate for consumer CD players and 48 kHz 
for professional digital audio. This, in turn, has given rise to sample-rate conversion ("SRC") systems for converting a 
stream of digital data at one sample rate to a stream of digital data at a different sample rate. However, the cost of exist- 
ing SRC systems is high. This imposes a severe constraint in designing more affordable digital audio products that 
15 apply to various source signals. 

[0003] Early implementations of SRC systems were done in a hybrid digital/analog domain. They were relatively 
simple, since all that is needed is a digital-to-analog (D/A) converter followed by an analog-to-digital (A/D) converter. 
The D/A converter runs at the input sample rate while the A/D converter is controlled by the output sample rate. If the 
output sample rate is lower, an analog anti-aliasing filter is provided between them. These three components are expen- 
se sive and consume a large amount of power, if designed for minimum signal degradation. 

[0004] Performing sample-rate conversion (SRC) in the digital domain has been a research/development topic for 
more than a decade. The article by R. E. Crochiere and L. R. Rabiner, "Interpolation and decimation of digital signals- 
A tutorial review," Proa IEEE, vol. 69, pp. 300-331 , March 1 981 , is an excellent reference for understanding fundamen- 
tal insights from early research results in this art area. Real-time, all-digital SRC systems are becoming more and more 
25 significant because digital processing of signals, such as voice, audio and video, appears to be increasingly dominant 
over traditional analog methods thanks to higher signal quality, rich features and the continually lowering cost of digital 
signal processing. 

[0005] Figure 1 shows a typical, all-digital SRC system 10 consisting of three basic building blocks: an interpolator 
(expander) 12, a high quality lowpass digital filter 14, and a decimator 16. The expander 12 takes an input stream of 
30 samples at one frequency, for example F s _ jn , and digitally produces a stream of digital samples at a higher rate that is 
an integer multiple, designated R in this example, of the input rate. Thus, the output of expander 12 is a stream of digital 
samples at a rate of 

F s _ out = R F s _ in . The decimator does the reverse. Thus, the decimator 16 takes an input stream of samples at one fre- 
quency, for example 

35 X = R F s _ jn , and digitally produces a stream of digital samples at a lower rate, divided by an integer division factor, des- 
ignated S in this example, of the input rate. This is referred to as decimation, or, alternatively, downsampling. Thus, the 
output of decimator 16 is a stream of digital samples at a rate of X/S = F s . out = (R/S)F s . jn . 

[0006] The reason for performing expansion followed by decimation is that the input sample rate and the output 
sample rate may not be a simple integer multiple of one another. The ratio of the interpolation factor R over the decima- 
40 tion factor S is the SRC ratio R/S where both R and S are positive integers. For R > S the SRC system is said to be 
operated in an SRC-UP mode, whereas an SRC-DOWN mode means R < S. 

[0007] R and S may be chosen to have large values in order to achieve higher quality SRC. However, when the val- 
ues of R and S are high, very-high order digital filters, usually finite impulse response ("FIR") digital filters, are neces- 
sary. Thus, as S and R increase, a greatly increasing amount of coefficient memory is required. For example, at a 
45 sample rate of 48 kHz an equi-ripple prototype FIR filter with transition bandwidth of Af = 4 kHz, passband ripples = 10* 
3 and stopband errors = 10~ 5 , has an order of approximately 128. The order of lowpass filter 14 in Figure 1 is 128R for 
an interpolation filter (R is the interpolation factor up to over 1000). 

[0008] It is actually quite difficult or even impossible as a practical matter to design an equi-ripple FIR filter with tens 
of thousands of orders. In fact, the major difficulty encountered in existing audio SRC system, implemented either in an 

so application-specific integrated circuit (ASIC) or on a programmable digital signal processor (DSP) such as a Texas 
Instruments TMS320 series DSP chip, seems to be large memory size and high computational complexity. For exam- 
ple, an ASIC described in an article by R. Adams and T. Kwan, "A stereo asynchoronous digital sample-rate converter 
for digital audio," IEEE J. Solid-State Circuits, vol. 29, pp.. 481-488, April 1994, needs tens of kilobytes of memory to 
store just a fraction of the nearly 10 million filter coefficients used. 

55 [0009] Moreover, filter coefficient interpolation, which is performed to generate thousands of sets of required poly- 
phase filter coefficients in real time, expends significant computational power which is provided by a hardware multiplier 
plus an accumulator. A similarly difficult situation is also encountered when a programmable DSP chip is employed. For 
example, see the article by S. Park et af., "A novel structure for. real-time digital sample-rate converters with finite pre- 
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cision error analysis," Proa Int. Conf. on Acoust, Speech and Signal Processing, pp. 3613-3616, Toronto/1991. Several 
kilobytes of memory are employed in these SRC systems for filter coefficients alone, in addition to their computational 
complexities falling in the neighborhood of 10 MIPs for one channel of high quality audio. 

[0010] Therefore, attempts have been made using window techniques, for example using a Kaiser window, to 

s design extremely high-order FIR filters. By using an interpolation technique in calculating required coefficients in real- 
time, a single-stage SRC filter system, such as the system 10 shown in Figure 1, has become closer to practical, and 
hardware implementation examples have been reported. Examples may be found, e.g., in U.S. Patent Nos. 4,780,892, 
4,564,918, 4,825,398 and 4,748,578. However, these implementations fall short of the desired efficiencies allowing 
their utilization in affordable digital audio products for consumers. 

10 [001 1 ] Other attempts at avoiding the use of high-order filters make use of special functions such as Lagrange pol- 
ynomials or B-spline functions. See, for example, T. O. Ramstad, "Digital methods for conversion between arbitrary 
sampling frequencies," IEEE Trans. Acoust, Speech and Signal Processing, vol. ASSP-32, pp. 577-591, June 1984, for 
an article on the former, and S. Cucchi et al., "DSP implementation of arbitrary sampling frequency conversion for high 
quality sound application," Proc. Int. Conf. on Acoust, Speech and Signal Processing, pp. 3609-3612, Toronto, 1991, 

is for the latter. These methods, however, all have the drawback of requiring a very large number of computations. 

[0012] It is known that multistage decimation or interpolation filters are generally more efficient than single-stage 
filters, in terms of computational complexity. It appears that the same conclusion also holds on memory requirements 
of multi-stage filters over single-stage versions. It would therefore be desirable to have a multi-stage scheme employing 
a far smaller memory than required in the prior art to store some of the SRC filter coefficients, and have an accompa- 

20 nying arrangement for efficiently calculating the rest of the filter coefficients, in real time. 

[001 3] Therefore, it is an object of this invention to provide an efficient multistage multi-rate filter. It is also an object 
of the present invention to provide a multistage SRC filter that is more efficient in both computational and memory 
requirements than prior art multistage SRC filter implementations. It is a further object of the present invention to pror 
vide a multistage SRC that represents a balance of resource considerations. 

25 '* 

Summary of the Invention 

[0014] In accordance with the present invention there is provided a method for providing a sample-rate conversion 
("SRC") filter on an input stream of sampled data provided at a first rate, to produce an output stream of data at a sec- 

30 ond rate different from the first rate. The input stream of sampled data is operated on with a first low-order interpolation 
filter routine to produce a first stream of intermediate data. The first stream of intermediate data is operated on with a 
first simplified interpolation filter routine, having a substantially small number of operations to calculate the coefficients 
thereof, to produce a second stream of intermediate data. The second stream of intermediate data is operated on with 
a first decimating filter routine to produce the output stream of data. 

35 [0015] Implementations of the present invention can exhibit excellent characteristics such as an extremely low 
memory requirement, simple design and implementing procedure and moderate computational complexity. 
[0016] Filter performance of the new implementations can easily reach 100 dB signal-to-noise-ratio (SNR) level. 
The simplicity inherent in the present invention allows hardware implementation, by using either an ASIC or a program- 
mable DSP chip, to be easier and more straightforward than heretofore. 

40 [0017] These and other features of the invention will be apparent to those skilled in the art from the following 
detailed description of the invention, taken together with the accompanying drawings. 

Brief Description of the Drawings 

45 [0018] 

Fig. 1 is a block diagram showing a typical, all-digital SRC system; 
Fig. 2 is a block diagram showing a preferred embodiment of the present invention; 
Fig. 3 is a graph showing the frequency response of a stage-3 filter, with r=3; 
so Fig. 4 is a more detailed diagram of a preferred embodiment of the present invention; 

Fig. 5 is a flow chart showing the general method for indexing used in the preferred embodiment of the present 
invention; 

Fig. 6 is a block diagram of a simplified model for explaining the indexing used in the preferred embodiment of the 
present invention; 

55 Fig. 7 is a block diagram of a more detailed model for explaining the indexing used in the preferred embodiment of 

the present invention; 

, Fig. 8 is a schematic diagram depicting a data output stream as a series of vectors; and 
Fig. 9 is a diagram of an SRC filter that utilizes al-A modulator. 
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Detailed Description of the Preferred Embodiment 

[001 9] The four stage SRC filter structure 20 depicted in Figure 2 is a block diagram of a preferred embodiment of 
the present invention. It represents an excellent balance between easy design procedure, reduced memory/computa- 
tion needs, and simple system control flow. The SRC filter structure 20 comprises a stage-1 expander/filter 22, a stage- 

2 expander/filter 24, a stage-3 expander/filter 26 and a stage-4 decimator 28. The stage-4 decimator 28 can be a simple 
factor of S prior art decimator, and so the discussion that follows focuses on how best to apply the principles of the 
present invention to design efficiently embodiments of the first three filters in stages 1 , 2 and 3, and how to design them 
for optimum performance. As an initial matter in that regard, only FIR filters are utilized in the embodiments disclosed 
herein, because they offer linear phase responses and efficient polyphase implementations in multirate processing. 
However, infinite impulse response ("IIR") digital filters may be used, if desired. 

[0020] Further, while the embodiment of the present invention disclosed in detail herein utilizes three expander/filter 
stages, this is only preferred, but not required. Two of such stages, or, alternatively, greater than three such stages may 
be used in the application of the principles of the present invention, and still remain within the scope of the present 
invention as defined by the appended claims. 

[0021 ] Now, returning to the preferred embodiment, the stage-1 expander/filter 22 performs an interpolation of fac- 
tor R-| . The value of R 1 is relatively small so that the transition bandwidth of the stage-1 expander/filter 22 is still reason- 
ably large. The reason for this is that a large value of , causing a high output sample rate at the stage-1 expander/filter 
22, can lead to an extremely narrow transition band that, in turn, proportionally results in a very high-order for stage-1 
filters, which is undesirable. It is desirable to maintain the order of the stage-1 filter between 100 and 800. 
[0022] With a small vale of R 1 , the output sample rate in the stage 2 expander/filter 24 is relatively high. However, 
even though the output sample rate is significantly higher in the stage-2 expander/filter 24 than in the previous stage! 
the transition band of the stage-2 expander/filter 24 can be more relaxed than that of stage 1 since a large number of 
the image bands have been already removed by the stage-1 expander/filter 22, preferably at least half of the image 
bands. The selection of R 2 is not arbitrary, however, because it is desirable to simplify the stage-3 expander/filter 26 
given that it is running at a rather high sample rate, R 1 R 2 F s . in . Thus, the choice of R 2 depends on how the stage-3 
expander/filter 26 is designed. 

[0023] The design aspects of the stage-3 expander/filter 26 occupy a significant role in the design of the preferred 
embodiment of the present invention. A key aspect is that the stage-3 expander/filter 26 is a relatively simple filter. Sim- 
ple-filters have trivial filter coefficients, and the simplest non-zero filter coefficients are simply ones. Further, the stage- 

3 expander/filter 26 is cascaded. As a general design matter, it should be appreciated that a trivial -coefficient-based fil- 
ter can perform acceptably if it is cascaded, even if only two times. 

[0024] The underlying principles for the stage-3 expander/filter 26 are now described. Let the stage-3 expander/fil- 
ter 26 be a cascade of r all-one-coefficients filters that are described by: 



[0025] Notice that, for simplicity, a normalization factor has been omitted in Equation (1). Using the filter described 
by Equation (1) in stage 3 substantially removes a significant hurdle that has existed heretofore in designing SRC sys- 
tems: huge storage requirements for filter coefficients. As is shown below, the coefficients of H 3 (z) can be evaluated 
based on very simple equations when r is small, say r=2 or 3. The magnitude response of H 3 (z) is easily evaluated as: 



[0026] The function of H 3 (z) is to eliminate the image bands introduced by padding (i.e., inserting) zero-value sam- 
ples between input samples of stage 3. Furthermore, the monotonic property of H 3 (z) between DC and the frequency 
<d 0 =2tt/R 3 , at which the filter frequency response reaches its first null value, allows the evaluation of Equation (2) at 
just two digital frequency points to determine the performance of the filter. 
[0027] The first point is 




Eq. (1) 




Eq. (2) 
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where f p denotes maximum frequency of base band signals, and R=R ,R 2 R 3 is the interpolation factor of the whole 
SRC system. At : 
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[0028] This expression gives the maximum droop in the passband of the stage-3 expander/filter 26. Thus 
is o> 2 =a> 0 -7t/R represents the edge frequency of the first stopband of the stage-3 filter, and the corresponding frequency 
response at g>2 is: 
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Eq. (4) 
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Note the following useful information that is evident from the above expression: 

Increasing the value of R 1 R 2 can greatly improve the filter performance of stage 3 in both passband and stopband. 
30 In other words, the first image band is located further away from base band as R 1 R 2 increases. 

For large values of R 3 , which indeed is the case in many SRC applications, the frequency responses at n<o t and 
qg>2 asymptotically become independent of R 3 . 

[0029] To quantitatively examine the filter performance of stage 3 at both G) 1 and a>2 several useful cases are eval- 
35 uated and shown in Table 1 with these parameters: R 3 = 256 (2 8 ), f p = 20 kHz and F s in = 48 kHz. 
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Stage-3 Filter Performance 


Ri 


R 2 


H(e^')(dB) 


H(ei^) (dB) 


2 


4 


-0.087a- 


-23.6r 


2 


8 


-0.022r 


-29.8r 


2 


16 


-0.006r 


-36.0r 


2 


32 


-0.001 r 


-42. ir 


2 


64 


-0.0003r 


^48. 1r 


2 


128 


-0.0001 r 


-54.2r 



[0030] A value of r=3 is an optimally balanced selection, and is considered preferred. It will be clear from an inspec- 
tion of Table 1 that a value of r=1 does not provide nearly adequate filter performance in most applications, although 
the filter coefficients are indeed trivial. 
55 [0031 ] With r=2, the coefficients of H 3 (z) still remain quite simple: 



h 3 (n) = n 1 < n > R 3 



Eq- (5) 
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Notice that the other half (R 3 - 1) of the filter coefficients (R 3 +1 < n > 2R 3 - 1) can be obtained from symmetry Although 
this filter provides good filter performance, over 1 00 dB stopband rejection for digital processing of audio signals is often 
desired. This requires the value of R 1 R 2 to be at least 256. It will be recalled that it is preferable to use R 1= 2 for the 
reason previously described. If R-j is selected to be 2, then to maintain the value of R 1 R 2 as at least 256 R 2 must be 
128 or even larger. Therefore, the sample rate at the output of the stage-2 expander/filter 24 reaches a' minimum of 
256F sJn . That is over 10 MHz when F sJn =44.1 kHz or higher. Also, for R 2 =128, the order of the stage-2 expander/filter 
24 may well increase to over a thousand in order to retain high filter performance. An order of one thousand means 
large amounts of memory are required, which is undesirable. It is desirable to maintain the order of the staae-2 filter 
between 100 and 1200. y 

[0032] With r=3, a high-quality filter is easily achieved even though R 1 R 2 has a moderate value. Figure 3 is a chart 
in which the horizontal axis represents frequency in Hertz and the vertical axis represents magnitude of response in dB 
showing such a filter's frequency response 30 within a frequency range including the first stopband 32 that rejects the 
first image band. The first stopband 32, or equivalent!* image band, starts at 1512 kHz and ends at 1560 kHz with the 
first null sitting at 1536 kHz. 

[0033] The filter coefficients when r=3 can be evaluated by the following expression; 

n(n + 1) 

2 J 1 < n > R 3 - 1 

Eq. (6) h 3 (n) = i 

^%^--(R 3 --nX2R3--^n) R 3 S n ^ ^L-l 



[0034] Notice that R 3 has been assumed to be an even number, and again only half of the symmetric coefficients 
are described in Equation (6). The end summing index in Equation (6) simply changes to (3R 3 -1)/2 if an odd integer of 
R 3 is chosen. The evaluating procedure to obtain those coefficients is sufficiently straightforward to allow computation 
in real time. Thus, each set of polyphase components of H 3 (z) consists of at most three integers that require a total of 
only three integer multiplications and several additions/shifts. The largest coefficient of H 3 (z) is given by 3R 3 R 3 /4=3x2 14 
for R 3 =2 , requiring merely 16-bit precision, which, in turn, requires only 16-bit registers. 

[0035] Given r=3, R 2 can be chosen as 8 or 1 6 for audio SRC applications. In the following design examples R 2 =1 6 
and R 3 =256. If strict passband performance is required a simple filter can be employed with several taps to pre-distort 
the base-band signals such that the droop in the passband introduced by the stage-3 expander/filter 26 is well compen- 
sated. 

[0036] it is well known that filters with fixed coefficients can be implemented more efficiently than those having pro- 
grammable coefficients. This is particularly true when implementing filters in ASICs. Therefore, when implementing 
H 3 (z), since the coefficients are determined by the value of R 3 and calculated in real time, the design must be general 
enough to accommodate every possible set of filter coefficients that will arise in a particular implementation In any 
event, it is desirable to maintain R 3 between 100 and 10,000. 

[0037] Upon having chosen the value of R 2 the stage-2 expander/filter 24 can easily be made with fixed coeffi- 
cients. Those fixed coefficients can further be represented by canonical signed digits (CSD) form to greatly simplify 
computations by using add/shift operations rather than normal multiplications. Even by using a programmable DSP 
chip, fixed coefficients are easier to handle than time-varying coefficients. 

[0038] In implementing an embodiment of the present invention particular attention should be made in designing 
the first stage filter. In this regard, it should be noted that making distinct the considerations between SRC-UP and SRC- 
DOWN modes has so far been ignored in the discussion herein. This will now be discussed. The care that should be 
exercised in designing efficient stage-1 filters that are able to handle both UP and DOWN modes highlights an important 
design consideration when seeking the high efficiency achievable by the inventive SRC design methodology disclosed 
herein. 

[0039] Some problems are encountered in design when UP modes are involved. They are discussed below The 
function of the stage-1 filter simply erases the image band caused by the first interpolator. Therefore the stage-1 
expander/filter 22 can have fixed coefficients. But in DOWN modes, the stage-1 expander/filter 22 must'remove high 
frequency components of the input signals in addition to cutting the image band off. 

[0040] Two embodiments of the present invention are presented herein for solving this problem, in the form of two 
methods. In the first method, several stage-1 filters may be provided, each of which is suitable to a specific DOWN 
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mode. The selection from those pre-set filters is easily done by feeding an index that points to a correct initial memory 
address of stored f ilter coefficients. If desired, those pre-computed coefficient sets can be stored in a host system such 
as a personal computer to save the memory on an ASIC or DSP chip. 

[0041 ] The second method calls for a real-time computing scheme, similar to H 3 (z), that generates filter coefficients 
5 based on a generic set of filter coefficients and still meets the filter specifications. Note that this computing procedure 
needs to be executed only once when the SRC ratio is changed, and stays on until a new SRC ratio is invoked. None- 
theless, the first of these two methods is considered preferred, and that method is utilized in the following design exam- 
ples. 

[0042] Designing filters that can be employed for converting sample rates between 44.1 and 48 kHz, based on the 
10 above design guidelines, will now be considered. 

A. Stage-1 Filters 

[0043] Stage-1 filters having an order of 1 26 provide adequate filter performance and require acceptable amounts 
is of memory and computational complexities. In the UP mode, i.e., the 44.1 > 48 kHz case, a 126th-order FIR filter can 
provide over 1 00 dB stopband rejection and maintain up to 20 kHz passband signals. It is also possible to employ a half- 
band FIR filter if a somewhat greater number of aliasing errors can be tolerated in the transition band from 20 kHz to 24 
kHz. In half-band FIR filters one of the two polyphase sub-filters is a delay element. Therefore half-band FIR filters need 
nearly 50% lower computations and coefficient memory in comparison to normal FIR filters with the same filter orders. 
20 To further reduce computations, symmetry in filter coefficients can be used. Indeed, by selecting an even-order, half- 
band FIR filters have polyphase components with symmetrical coefficients. 

[0044] In the DOWN mode, however, a 126th-order FIR filter can only maintain a passband of about 18.5 kHz if the 
100 dB stopband rejection requirement applies. 

[0045] Actual filter coefficients for a Stage-1 filter used in a Matlab simulation of an embodiment of the present 
25 invention are found in Appendix A. The designer may wish to select other coefficients to accommodate specific appli- 
cation constraints and/or performance objectives. Optimal quantization of filter coefficients may be considered by the 
designer, as well. All of such variations are well within the purview of those of ordinary skill in this art area. 

B. Stage-2 Filters 

30 

[0046] Only the stage-2 expander/filter 24 should have fixed coefficients that are independent of the SRC mode. To 
be consistent with the filter performance in the stage-1 expander/filter 22 configured as described above, a 143rd-order 
FIR filter should be used. The order number 1 43 offers an advantage in that every polyphase component has the same 
number of coefficients, that is, nine, when R 2 =16. Although the polyphase components may no longer have symmetric 

35 coefficients, half of the polyphase components have coefficients that are mirror-symmetric with their counterparts in the 
other half of the polyphase components. Such a property may be exploited to save coefficient memories. 
[0047] Actual filter coefficients for a Stage-2 filter used in a Matlab simulation of an embodiment of the present 
invention are found in Appendix B. The designer may wish to select other coefficients to accommodate specific appli- 
cation constraints and/or performance objectives. Optimal quantization of filter coefficients may be considered by the 

40 designer, as well. All of such variations are well within the purview of those of ordinary skill in this art area. 

[0048] In this connection, the designer may wish to consider for stage-2 filters a cascade of several half -band FIR 
filters having half zero filter coefficients, which results in a need for only half of the computations otherwise, as well as 
half of the memory requirement. The designer should keep in mind, however, that a drawback of using a cascade of 
several half-band FIR filters is the need for a more complicated indexing scheme than that described in the following 

45 sections. 

C. Stage-3 Filters 

[0049] This filter has already been discussed in detail in the previous section. Note, however, that it might be pre- 
so ferred to make this filter fully programmably controlled by selecting R 3 in real time. Remember that larger values of R 3 
require higher dynamic ranges in filtering computations. 

[0050] Since the impulse response of this filter is very smooth, it is possible to compress the coefficients. Piecewise 
linearization may be used, for example. However, since the coefficients are already easily calculated, and since the fil- 
tering operation does not use neighboring coefficients in the computation of any given output sample, it is not clear that 
55 significant savings would result from compression, but the implementation enhancement is suggested in the interest of 
full disclosure. 
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D. Implementation Architecture 

[0051 ] A more detailed diagram of a preferred implementation of the new SRC filter is shown in Figure 4. The SRC 
filter 40 includes three expander/filter stages 42, 44, 46, and a decimator stage 48. The stage-1 expander/filter 42 
includes two polyphase sub-filters 50, 52, corresponding to R 1= 2, each receiving the input x(n) provided at a sample 
rate of F sJn . The outputs of polyphase sub-filters 50, 52, are selected at a rate of 2F sJn , as shown by switch 54, and 
provided at that rate as the input x 1 (n) to the stage-2 expander/filter 44. 

[0052] The stage-2 expander/filter 44 includes sixteen polyphase sub-filters 60, 62, 64, ... 66, corresponding to 
R 2 =16, each receiving the input x t (n). The outputs of polyphase sub-filters 60, 62, 64, ... 66, are selected sequentially 
and cyclically at a rate of 32F sJn , as shown by switch 68, and provided at that rate as the input q(n) to the stage-3 
expander/filter 46. 

[0053] The stage-3 expander/filter 46 includes R 3 polyphase sub-filters 70, 72, 74, ... 76, each receiving the input 
q(n). The outputs of polyphase sub-filters 70, 72, 74, ... 76, are selected sequentially and cyclically at a rate of 
32R 3 F sJn , as shown by switch 78, and provided at that rate as the input to the stage-4 decimator 48. The output of dec- 
imator 48 is the output y(n) of the SRC filter 40. 

[0054] In practice, the stage-4 factor-S decimating must be merged into the three previous stages to avoid any 
redundant computations that are not related to output samples at a required sample rate. This is accomplished in the 
preferred embodiments herein by a novel indexing technique, explained below. The resulting memory size for the SRC 
filter coefficients is only 136 words, which is more than an order of magnitude reduction in comparison to the prior art 
SRC systems described in the article by R. Adams and T. Kwan, and the article by S. Park era/., cited above. On the 
other hand, the number of multiplications is 97 per output sample, a moderate quantity for a high quality SRC system. 
[0055] An SRC procedure for converting samples provided at a rate of 48 kHz to a rate of 44.1 kHz will now be 
explained in detail, in conjunction with a Matlab listing implementing such procedure. This procedure includes a novel 
indexing scheme alluded to above. It will be recalled that it was pointed out that major increases in storage and process- 
ing efficiencies are obtained in the preferred embodiments by way of such novel indexing scheme. This indexing is best 
understood by working backwards from the output sample. Figure 5 is a flow chart showing the steps in determining the 
current output sample. Thus, the first step 80 is to determine which stage-3 output sample will be the current SRC out- 
put sample. The second step 82 is to determine which stage-2 output samples will be used to form the appropriate 
stage-3 output sample. The third step 84 is to determine the coefficients of the stage-3 filters that will be applied to the 
stage-2 output samples to form the needed stage-3 output sample, and to compute those coefficients. The fourth step 
86 is to determine the stage-1 output samples needed to form the needed stage 2 output samples. The fifth step 88 is 
to determine the stage-2 filter coefficients needed to form the needed stage-2 output samples, and then to retrieve 
those coefficients from storage. The sixth step 90 is to determine the input samples needed to form the needed stage- 
1 output samples. The seventh step 92 is to determine the stage-1 filter coefficients needed to form the needed stage- 
1 output samples, and then to retrieve those coefficients from storage. The eighth step 94 is to perform the necessary 
filtering operations for stages 1 , 2, then 3. The ninth step 96 is to output the sample and return to the first step 80. 
[0056] The process described above generally, in connection with Figure 5 will now be explained in more detail. 
First, Figure 6 is a simplified block diagram showing an SRC filter 140, based on the SRC filter 40 of Figure 4, for the 
case where R 3 = 147, S = 5120 (this is for the case of 48 kHz to 44.1 kHz conversion). The input sample rate F s in is 
referred to here, and hereinafter in this document, as F s , to simplify notation. Thus, the output of the stage-1 
expander/filter 142 is provided at a rate of 2F S . The output of the stage-2 expander/filter 144 is provided at a rate of 
32F S . The output of the stage-3 expander/filter 146 is provided at a rate of 147*32F S . Decimation is performed by 
selecting output samples from the output of the stage-3 expander/filter 1 46 at the rate of 1 60*32F S , as shown by switch 
148. The resulting output y(n) is a stream of samples provided at a rate of (HZ )F S . Since out of 147*32 = 4704 samples 
on, y (im ) ~ 1 s a m Ples are actually provided as output samples, many samples output from the stage-3 expander/filter 
1 46 are discarded. The preferred embodiment takes advantage of this and only calculates those samples actually pro- 
vided as output samples. Some additional calculation is performed in order to accomplish this, but since the reduction 
in sample calculations is by a factor of over 4,000, this additional calculation results in a significant reduction in memory 
and processing requirements. 

[0057] Figure 7 is a block diagram of the SRC filter 140 of Figure 6 expanded in detail to assist in understanding 
the selection process back through the stages, as explained above in connection with Figure 5. 
In implementing the functionality of Figure 7, Figure 8 will also be referred to, to show how to index the SRC system 
such that only necessary computations are carried out. 

[0058] Note first that if the indexing scheme were not being used, a data stream at 32*147*Fs would be seen prior 
to the downsampling, or decimation, step. However, the downsampling step provided in the preferred embodiment of 
the present invention takes only one out of every 32*160 of these samples. Referring to Figure 8, consider the output 

of stage 3 prior to downsampling to be vectors 1 50, 152 1 54, 1 56, of 1 47 samples coming at a rate of 32Fs. This is 

true since there are 147 polyphase filters in stage 3 146 (Fig. 7). To accomplish the downsampling, it is determined 
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which outputs, at which time steps, actually become outputs. After such determination, the method of the referred 
embodiment works backwards, as it were, to find out which samples must be used with which filter coefficients in the 
various stages to compute this output. 

[0059] Table 2, below, shows which samples are used if there is a 147*32*Fs Hz data stream being downsampled 
5 by 160*32. 



Table 2 



Output 


Time 


Filter 


1 


1 


1 


32*160+1 


35 


123 


2*32*160+1 


70 


98 


3*32*160+1 


105 


73 


4*32*160+1 


140 


48 



This data are provided by determining which filter outputs provide desired samples, at which time steps For example, 
20 the 4 th output sample is sample 3*32*160+1 prior to downsampling. This is the output of filter 73 at time step 105. 
Notice that, since the output of stage 2 is running at 32Fs, the output of the stage 2 that is occurring at the 105 th time 
step is exactly that 105 th sample. 

[0060] Table 3, below, shows the general equation used to compute the output samples needed, the stage 2 sam- 
ples that correspond to those outputs, and the stage 3 filters that will be involved in computing the final output. 

25 

Table 3 



Output 


Stage 2 
Sample 


Stage 3 
Filter 


1 


1 


1 


32*160+1 


35 


123 


2*32*160+1 


70 


98 


3*32*160+1 


105 


73 


4*32*160+1 


140 


48 


0=(i- 

1)*32*160+1 ! 


TxWo/14 71 


Fi=0- (T- 
1) *147 



[0061] Table 4 ( below, shows the stage 2 filters and stage 1 samples needed to generate the downsampled data 
stream. 



55 
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Table 4 



Stage 2 
Sample 


Stage 1 
Sample 


Stage 2 
Filter 


1 


1 


1 


2 


1 


2 








16 


1 


16 


17 


2 


1 








m 


T c =U/16l 


F 2 =m-(T- 
1)*16 



[0062] The following Matlab code shows the specific case of 48kHz to 44.1kHz conversion. The efficient indexing 
scheme is applied only to the 2 nd and 3 rd stages, but could easily be extended to the 1 st stage as well And while this 
Matlab code shows only the case of 48kHz to 44.1 kHz conversion, the teachings of this section can easily be extended 
to indexing schemes for any conversion rate. 
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Matlab Listing: 

^Assumes 48kHz in, 44.1kHz out. 

%1 Channel - multiple channels are just replicas 
function src 

%stage 1 does not implement efficient indexing 
%Zero pad 

ss = zeros (1, 2*max (size (in) ) ) ; 
ss (1:2: end) = in; 

%(for "in", start at 1, skip by two, continue to the end) 
%Filter with stage 1 filter 

ss = 2*f liter (H1_B, H1_A, ss) ; . ■ • 

%(the H1JB and Hl_A vectors are read from storage) 

%make polyphase filters, i.e., generate the FIR coefficients 
for stage 2 

%from non-polyphase stage 2 filter 
for i = 1:16 

H2(i, :) = 16*H2 (i : 16:end) ; 

end 

%Compuation for stages 2 and 3 
for i = 1 : f ix (147*numpts/160) -9 

%subtract 9 since the first nine input samples are discarded 

%to simplify initialization. 
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%stage 3 indicies 

O = (i - 1 ) +160*32 + 1; 
T2 = ceil (0/147); 
F2 = O - (T2 - 1) + 147; 

Istage 3 coefficients 

%function call; go to function "A", below 
C = coef (F2) ; 

%compute necessary stage 2 data 
for j = 1:3 

%stage 2 indicies 

m = n + j - 1; 

Tl = ceil (m/16) ; 

Fl = m - (Tl - 1) +16. 

Tl = Tl + 9; 

%adding 9 insures that the index will not go beyond the 
%file beginning, this is equivalent to discarding the 
Sfirst 9 input samples. 

^perform stage 2 filters, i.e., convolve with stage 

output 

vec(j) = 147*B2(F1, :) + ss (Tl : -1 : Tl-8 ) ■ ; 

end 

^perform stage 3 filter 
in (i) = C+vec f ; 

end 



out(l:fix(147+numpts/160)-9) = in (1 r f ix (147+numpts/160) -9) ; 
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20 



30 



35 



^Compute Stage 3 Filter Coefficients 
% function "A" : 
function C = coef (m) 

C = zeros (1/ 3) ;' 



%calculate Eq. (6) (from text, above) 
for i = 1:3 
15 n - (i - 1) * 147 + m; 

if n >= 1 & n < 147 

C(i) = n* (n+l)/2; ' 
elseif n >= 147 & n < 2*147 

C(i) = 147* (147+1) /2 + (n - 147) * (2*147 - 1 - n) ; 
elseif. n >= 2*147 & n < 3*147-1 
25 c(i) = (3*147 - n - 1)* (3*147 - n)/2; 

end 



end 

Iprepare for convolving 
C = fliplr (C) /3176523; 



40 VI. £-A OVERSAMPLING ANALOG-TO-DIGITAL AND DIGITAL-TO- ANALOG CONVERTERS WITH EMBEDDED 
SRC FUNCTIONS 

[0063] A further implementation variation will now be discussed, that is, the use of the interpolation filter in a £ - A 
oversampling D/A converter in the place of the stage- 1 expander/filter 22 and the stage-2 expander/filter 24 of the SRC 

45 system 20 of Figure 2. In this way, an SRC function may readily be incorporated into a Z - A oversampling D/A converter. 
[0064] Note, however, that whereas a normal SRC filter does not have oversampled output values, a £ - A oversam- 
pling D/A converter requires oversampled signals to feed the digital modulator. Therefore, it is important not to simply 
add a £ - A D/A interpolation filter to an SRC construded in accordance with the principles of the present invention, 
because, for example, the stage-3 expander/filter 26 of Figure 2 may have to generate many more output samples in 

so such a case. Inventive prindples to accomplish such a combination will now be described. 

[0065] First, it is advantageous to employ a modified SRC filter that produces an output sampled at a sample rate 
eight times its Nyquist sample rate. This takes advantage of a known interpolation filter technique in a £ - A oversam- 
pling D/A converter, that is, to upsample signals to an eightfold sample rate using a lowpass filter followed by a sample- 
and-hold drcuit that boosts the sample rate to a higher level, for example to 128 times the original sample rate. 

55 [0066] Implementing such an embedded SRC filter adds considerable computational complexity. To reduce the 
computational complexity, two optional implementations may be employed. First, note that the interpolation filter per- 
formance in a S - A oversampling D/A converter can be lower than that of a normal SRC filter, since high-frequency 
shaped noise is introduced by a digital modulator. Therefore, a filter with an 80-dB stopband rejection may be employed. 
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As a result, for a selection of R 2 =8 the performance of H 3 (z) is now higher than adequate. 

[0067] Figure 9 is a diagram of a decimation filter 1 62 that follows a E - A modulator 1 60. The modulator 1 60 works 
at a rate Of F mod =3.072 MHz that is typically an oversampled rate of 64 to 384 times the first sample rate F S1 . A fourth- 
order SINC filter 164 is employed to downsample the modulator output signals to a rate of 16F S1 . Then, a simple filter 
1 66 such as the stage-3 expander/filter 26 of Figure 2, with an order of N=r(S-1) , is used to convert the sample rate at 
the SINC filter 164 output into a rate of 16 times the final sample rate F S2 . 

[0068] This simple filter 1 66 performs an SRC function with a ratio of R/(S=160) . It should be pointed out that the 
order of the filter is based on the decimation factor S rather than the interpolation factor R. The reason for the modifica- 
tion lies in the fact that there is now an anti-aliasing requirement, while such a condition is not required in other embod- 
iments described above. In this embodiment, however, R has been optimized to have only three possible values: 144, 
147 or 160. Note that when R=1 60 a bypass function results. The final stage 168 is a 16:1 decimator. 
[0069] All possible SRC cases are summarized in Table 5. From an inspection of Table 5 it is apparent that a wide 
range of frequently used audio sample rates may be converted using this approach. 



Table 5 



Sample Rates in A £ - A A/D with SRC 


M-| 


16F S1 


R/(S=160) 


16F S2 


M 2 


F S 2 


4 


768k 


1 


768k 


16 


48k 


4 


768k 


147/160 


705.6k 


16 


44.1k 


6 


512k 


1 


512k 


16 


32k 


8 


384k 


147/160 


352.8k 


16 


22.05k 


10 


307.2k 


1 


307.2k 


16 


19.2k 


12 


256k 


1 


256k 


16 


16k 


12 


256k 


144/160 


230.4k 


16 


14.4k 


16 


192k 


147/160 


176.4k 


16 


11.025k 


20 


153.6k 


1 


153.6k 


16 


9.6k 


24 


128k 


1 


128k 


16 


8k 


24 


128k 


144/160 


115.2k 


16 


7.2k 



[0070] The design example shows that an inband SNR of over 90 dB can easily be achieved by using the decima- 
tion/SRC filter in Figure 9 with four typical half-band FIR filters performing a decimation of M 2 =16, just like those filters 
employed in normal X - A oversampling A/D converters. 

[0071] A compact filter procedure is now provided that implements the 'stage-3' filter H 3 (z) 166 in Figure 9 as fol- 
lows: 

// The 'Stage-3' Filter Calculation Program 

// The SRC ratio is R/S: R is the interpolation and S the 
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decimation ratio. 

// L=r(S-l)+l is the 'stage-3' filter's impulse response 
length. 

// M=Floor(L/R) is the maximum number of terms calculated 
for each output* 

// INPUT is an array of delayed input values. 
// INPUT [0 ] is the most recent input. 

// The RESAMPLER_INDEX register keeps track of when an 
output may be 

// calculated as well as providing the convolutional offset 
between the // input sequence and the filter impulse response. 
Filter coefficients can be 

// evaluated based on (5) or (6) with R 3 being replaced by 
S=160. 

1) Initialization 

RE SAMPLER^ INDEX = R; 
INPUT [all] = 0; 

2) do forever { 

wait for new input; 

RESAMPLER_INDEX = RESAMPLER_INDEX-R; 
while RESAMPLER_INDEX < R do { 
SUM =0; 

COEFFICIENT_INDEX = RESAMPLER_ INDEX + 1; 
for DAT A_ INDEX = 0 to M-l do ( 

SUM = SUM + INPUT [DATA_INDEX] 

FILT (COEFFICIENT_INDEX) ; 

COEFFICIENT_INDEX = COEFFICIENT_INDEX + R; 

if COEFFICIENT_INDEX > L then exit the inner 

loop; } 

y(n) = SUM; // y(n) is the output 
RESAMPLER_INDEX = RE SAMPLE R_JENDEX + S; } 

} 

[0072] Although the present invention and its advantages have been described in detail, it should be understood 
that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of 
the invention as defined by the appended claims. 

[0073] For example, a modification of the sample-rate conversion ("SRC") filter according to the invention receiving 
an SRC input and outputting an SRC output can be such that it comprises a first, low-order interpolation filter receiving 
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said SRC input; a second, low-order interpolation filter receiving the output of said first interpolation filter, said first inter- 
polation filter and said second interpolation filter being cascaded; a third, interpolation filter receiving the output of said 
second interpolation filter, and having a substantially higher order than said first interpolation filter and said second 
interpolation filter, and having a substantially small number of operations to calculate the coefficients thereof; a first dec- 
imating filter receiving the output of said third interpolation filter; wherein an optimized indexing procedure that automat- 
ically determines a subset of operations required to determine said SRC filter output is used to reduce computations 
involved in implementing said filters. 

[0074] Moreover the third interpolation filter advantageously can be such that it includes no normalization factor. 
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Appendix A 

Stage-1 Filter Coefficients 
-0 . 0 000 08 000080.31 
0.00002524520299 
0.00012749245978 
0.00023168895216 
0. 00018596926860 
-0.00004173798314 
-0.00022236139097 
15 -0.00009395165195 

0.00022798310620 
0.00028236300811 
-0.00011124263531 
-0.00045226922162 
-0.00015723875420 
0.00049028Q87929 
25 0. 000527 66243454 

-0.00029280537329 
-0.00086818706856 
-0.00017897003432 

30 

0. 00099041381460 
0. 00084653891504 

-0.00071764271163 
35 -0.00149994275425 

-0.00002733586373 
0 . 00183787295074 
0. 00114604821292 

40 

-0 . 00156196018316 
-0.00233378596337 
0.00050741782361 
45 0.00312720489457 

0.00123561352837 
-0.00303446960299 
-0. 00326181826590 
0.00172250243655 
0.00489531743842 



50 



55 
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0. 00079449580461 
-0.00535468499783 
-0.00403846539391 
0.00400977702922 
0.00708514672038 
-0.00066361611960 
-0. 00874702117957 
-0.00423331569757 
0.00790147137539 
0.00952642115799 
-0. 00388972863157 
-0.01351188782368 
-0.00312368468930 
0.01429804952811 
0.01194602726347 
-0.01029924266330 
-0.02038200385225 
0.00073482566492 
0. 02548598834348 
0. 01401306932399 
-0.02391373911549 
•0 . 03230583189130 
0.01201952643100 
0. 05147573458066 
0. 01540870560244 
0.06834265722374 
0. 07388677535303 
0.07990266018034 
0. 30676295531302 
0. 41556780905067 
0. 30676295531302 
0. 07990266018034 
0. 07388677535303 
0. 06834265722374 
0 . 01540870560244 
0 . 05147573458066 
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20 



0 .01201952643100 
-0.03230583189130 
-0.02391373911549 
0.01401306932399 
0.02548598834348 
0 .00073482566492 
-0.02038200385225 
-0.01029924266330 
0.01194602726347 

15 0.01429804952811 

-0.00312368468930 
-0.01351188782368 
-0.00388972863157 
0.00952642115799 
0.00790147137539 
-0.00423331569757 
25 -0.00874702117957 

-0.00066361611960 
0.00708514672038 
0.00400977702922 

30 

-0.00403846539391 
-0.00535468499783 
0. 00079449580461 

35 0. 00489531743842 

0.00172250243655 
-0.00326181826590 
-0. 00303446960299 

40 

0. 00123561352837 
■ 0..00312720489457 
0.00050741782361 

45 -0.00233378596337 

-0.00156196018316 
0.00114604821292 
0.00183787295074 
-0.00002733586373 
-0.00149994275425 



50 



55 
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-0 . 00071764271163 
0 . 00084653891504 
0 . 00099041381460 
-0.00017897003432 
-0. 00086818706856 
-0. 00029280537329 
0.00052766243454 
0.00049028087929 
-0. 00015723875420 
0. 00045226922162 
0.00011124263531 
0.00028236300811 
0. 00022798310620 
0. 00009395165195 
0. 00022236139097 
0. 00004173798314 
0. 00018596926860 
0. 00023168895216 
0. 00012749245978 
0 . 00002524520299 
0 . 00000800008031 
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Appendix B 

Stage-2 Filter Coefficients -0.00000656642335 
-0.00000829347347 
-0.00001299396513 
-0. 00001889619486 
-0.00002590574816 
-0.00003377059844 
-0.00004206107605 
-0.00005010239889 
15 -0.00005699368309 

-0.00006157598622 
-0.00006245406490 
-0.00005801156059 
-0.00004647574846 
-0.00002598701520 
0.00000529261223 
25 0. 00004 90 6099173 

0.00010672316139 
0.00017923064540 
0.00026690249991 

30 

0 .00036925440324 
0.00048483111855 
0.00061106542586 
35 0.00074416388891 

0.00087904339359 
0.00100932243685 
0.00112738671949 

40 

0.00122452909280 
0. 00129117601342 
0.00131719649461 
45 0.00129229867475 

0.00120650256525 
0.00105068232794 
0.00081715607669 
0.00050030613397 
0.00009720076660 



50 



55 
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-0 • 00039180669783 
-0. -00096253565435 
-0 . 00160652184989 
-0 .00231068193092 
-0.00305712544004 
-0.00382314240821 
-0.00458138407311 
-0.00530025283113 
-0.00594450502127 
-0. 00647606554720 
-0.00685503945042 
-0.00704089968951 
-0.00699381729954 
-0.00667609485928 
-0.00605365327878 
-0.00509751898954 
-0.00378525193285 
-0.00210225610565 
-0.00004291296422 
0 . 00238851614791 
0.00517726864880 
0 . 00829777990369 
0.01171373342216 
0. 01537847794018 
0.01923581503575 
0. 02322114491056 
0. 02726294492186 
0.03128454006421 
0.03520611288087 
0. 03894688805497 
0.04242741926498 
0.04557189914336 
0.04831041122226 
0. 05058104253984 
0.05233178003375 
0.05352212044397 
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0.05412433393826 

0.05412433393826 

0. 05352212044397 

0.05233178003375 

0.05058104253984 

0.04831041122226 

0.04557189914336 

0.04242741926498 

0. 03894688805497 

0.03520611288087 

0.03128454006421 

0. 02726294492186 

0.02322114491056 

0.01923581503575 

0. 01537847794018 

0.01171373342216 

0. 00829777990369 

0.00517726864880 

0.00238851614791 

0.00004291296422 

0.00210225610565 

0.00378525193285 

0 . 00509751898954 

0.00605365327878 

0.00667609485928 

0.00699381729954 

0.00704089968951 

0.00685503945042 

0.00647606554720 

0.00594450502127 

0. 00530025283113 

0.00458138407311 

0.00382314240821 

0.00305712544004 

0.00231068193092 

0.00160652184989 
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-0. 00096253565435 
-0 . 00039180669783 
0 . 00009720076660 
0 . 00050030613397 
0,00081715607669 
0.00105068232794 
0.00120650256525 
0.00129229867475 
0.00131719649461 
0 .00129117601342 
0 . 00122452909280 
0 . 00112738671949 
0,00100932243685 
0. 00087904339359 
0.00074416388891 
0.00061106542586 
0.00048483111855 
0. 00036925440324 
0.00026690249991 
0. 00017923064540 
0.00010672316139 
0. 00004906099173 
0.00000529261223 
-0.00002598701520 
-0.00004647574846 
-0.00005801156059 
-0. 00006245406490 
-0.00006157598622 
-0. 00005699368309 
0. 00005010239889 
0.00004206107605 
0.00003377059844 
0 . 00002590574816 
0 . 00001889619486 
0.00001299396513 
0.00000829347347 
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Claims 

1. A method for providing a sample-rate conversion ("SRC") filter on an input stream of sampled data provided at a 
first rate, to produce an output stream of data at a second rate different from said first rate, comprising the steps of: 

operating on said input stream of sampled data with a first low-order interpolation filter routine to produce a first 
stream of intermediate data; 

operating on said first stream of intermediate data with a first simplified interpolation filter routine, having a sub- 
stantially small number of operations to calculate the coefficients thereof, to produce a second stream of inter- 
mediate data; and 

operating on said second stream of intermediate data with a first decimating filter routine to produce said out- 
put stream of data. 

2. A method as in Claim 1 , further comprising the step of 

operating on said first stream of intermediate data with a second low-order interpolation filter routine to pro- 
duce a third stream of intermediate data; 

and wherein said step of operating on said first stream of intermediate data with a first simplified interpolation 
filter routine is performed by operating on said third stream, rather than said first stream of intermediate data. 

3. A method as in Claim 1 , further comprising the step of 

applying an optimizing indexing procedure in performing instructions of said routines so as to not execute 
instructions that do not generate intermediate data on which said output stream of data is based. 

4. A method as in Claim 2, wherein said step of operating with said first low-order interpolation filter routine is per- 
formed by interpolating with an FIR digital interpolation filter having an order of between 100 and 800, and that pro- 
vides said first stream of intermediate data at a rate R-j times said first rate, where R-j is an integer less than six. 

5. A method as in Claim 2, wherein said step of operating with said second low-order interpolation filter routine is per- 
formed by interpolating with an FIR digital interpolation filter having an order of between 100 and 1200, and that 
provides said second stream of intermediate data at a rate R 2 times said first rate, where R 2 is an integer between 
ten and thirty-six. 

6. A method as in claim 1, wherein said step of operating with said first simplified interpolation filter routine is per- 
formed by interpolating with a cascaded FIR digital interpolation filter having r filters in cascade relation, where r is 
an integer between 1 and 5, having coefficients that are all ones, having an order of R 3 -1 , and that provides said 
second stream of intermediate data at a rate R 3 times said first rate, where R 3 is an integer between 100 and 
10,000. 

7. A method as in Claim 2, wherein said step of operating with said first simplified interpolation filter routine is per- 
formed by interpolating with a cascaded FIR digital interpolation filter having r filters in cascade relation, where r is 
an integer between 1 and 5, having coefficients that are all ones, having an order of R 3 -1, and that provides said 
second stream of intermediate data at a rate R 3 times said first rate, where R 3 is an. integer between 100 and 
10,000. 

8. A method as in Claim 6, wherein said step of operating with said first simplified interpolation filter routine is per- 
formed by applying said second stream of intermediate data to a digital FIR filter according to the transfer function 
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H 3 (z)=f£ 1 z- k y 



A method as in Claim 7, wherein said step of operating with said first simplified interpolation filter routine is per- 
formed by applying said second stream of intermediate data to a digital FIR filter according to the transfer function 



H 3 (z)=[ Xz- k ] 

^k=0 ^ 



10. 



A method as in Claim 7, further comprising the step of applying an optimizing indexing procedure in performing 
instructions of said routines so as to not execute instructions that do not generate intermediate data on which said 
output stream of data is based. 



1 1 . A method as in claim 1 , wherein said second rate is R/S times said first rate. R being an interpolation factor and S 
being an decimation factor, said first low-order interpolation filter routine being a X-A modulator filter routine produc- 
ing said first stream of intermediate data at an oversampled rate between 64 and 384 times said first rate; 

said first simplified interpolation filter routine being a SINC filter routine having a fourth order to produce said 
second steam of intermediate data downsampled to a rate of sixteen times said first rate; 
modifying said second stream of intermediate data with a simplified interpolation filter routine having a sub- 
stantially small number of operations to calculate the coefficients thereof, by applying a cascaded FIR digital 
interpolation filter having r filters in cascade relation, where r is an integer between one and five, having coef- 
ficients that are all once, having an order of r(S-1), and that provides a modified second stream of intermediate 
data at a rate sixteen times said final rate; 

operating on said modified second stream of intermediate data with a decimation filter routine to produce said 
output stream of data. 

12. A sample-rate conversion ("SRC") filter receiving an SRC input and outputting an SRC output, comprising: 

a first, low-order interpolation (22) filter receiving said SRC input; 

a second, interpolation filter (26) receiving the output of said first filter, and having a substantially higher order 
than said first filter (22), and having a substantially small number of operations to calculate the coefficients 
thereof; 

a third, decimating filter (28) receiving the output of said second filter. 

13. A sample-rate conversion ("SRC") filter as in Claim 12, further comprising a fourth, low-order interpolation filter (24) 
disposed intermediate said first filter (22) and said second filter (26), receiving the output of said first filter (22) 
wherein said second filter (26) receives the output of said fourth filter (24), rather than the output of said first filter 

14. A sample-rate conversion ("SRC") filter as in Claim 12, further comprising means for automatically determining a 
subset of operations required to determine said SRC filter output, and for causing only said subset to be performed, 
such that the computations involved in implementing said filters are reduced in number. 



15, 



A sample-rate conversion ("SRC") filter receiving an SRC input and outputting an SRC output, comprising: 
a first, low-order interpolation filter receiving said SRC input; 

a second, low-order interpolation filter receiving the output of said first interpolation filter, said first interpolation 
filter and said second interpolation filter being cascaded; 

a third, interpolation filter receiving the output of said second interpolation filter, and having a substantially 
higher order than said first interpolation filter and said second interpolation filter, and having a substantially 
small number of operations to calculate the coefficients thereof; 
a first decimating filter receiving the output of said third interpolation filter; 

wherein an optimized indexing procedure that automatically determines a subset of operations required to 



26 



EP 1 039 636 A2 

determine said SRC filter output is used to reduce computations involved in implementing said filters. 

16. A sample-rate conversion ("SRC") filter as in Claim 15, wherein said third interpolation filter includes no normaliza- 
tion factor. 
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